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CO ' The Faddeev technique is employed to study the influence of both particle-particle and particle- 

hole phonons on the one-hole spectral function of le O. The formalism includes the effects of nuclear 
fragmentation and accounts for collective excitations at a random phase approximation (RPA) 
level. The latter are subsequently summed to all orders by the Faddeev equations to obtain the 
nucleon self-energy. The results indicate that the characteristics of hole fragmentation are related 
to the low-lying states of le O and an improvement of the description of this spectrum is required 
to understand the details of the experimental strength distribution. A first calculation to improve 
on this includes the mixing of two-phonon configuration with particle-hole ones. It is seen that 
this allows to understand the formation of excited states which cannot be described by the RPA 
approach. 
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I. INTRODUCTION 



The studies of spectral distribution in nuclei during the last decade has made available the best 
information to date regarding correlations in nuclear systems [jj. Experimentally, (e,e'p) reactions 
have provided precise information for spectroscopic factors of states close to the Fermi energy 0, 01 
in finite nuclei. Moreover, information on the tail of the spectral distribution at high missing energies 
and momenta, generated by short-range correlations (SRC), has recently become available [4]. The 
observed spectroscopic factors have been reproduced for selected nuclei, such as 48 Ca and 90 Zr j5|, and 
a particularly successful comparison between theory and experiment has been obtained for 7 Li 0. 
In all these calculations, the effects of configuration mixing at low energy was taken properly into 
account. Besides, the effects of SRC have been investigated for different systems and are known 
to give a constant reduction of the spectroscopic strength of 10-15%, for the fragments close to the 
Fermi energy 0, H, 0] • 

For the case of 16 O, there still exists a substantial disagreement between experiment and 
theory 0, 0, E3, EH The latter has been able to explain only a part of the quenching of 

spectroscopic factors for proton knockout from the p-shell. The results of Refs. [HI El show that 
the reasons for this discrepancy have to be looked for in terms of long-range correlations (LRC) and 
in particular in the couplings of single-particle (sp) motion to low-ener gy collective excitations. 



In this paper we report about the work done along this line in Refs. [12|, |13[ in order to approach 



* Proceedings of the "10 th international conference on nuclear reaction mechanisms", Varenna, 9-13 June 2003. 

^Electronic address: barbieri@triumf.ca URL: http://www.triumf.ca/people/barbieri 

t Electronic address: wimd@wuphys.wustl.edu URL: http://www.physics.wustl.edu/~wimd 



2 











\ 








o ♦ 


R (2p1h) 




+ 




R (2h1p) 















JLvaaaaa/^ 



FIG. 1: Diagrams contributing to the irreducible self-energy S*. The double lines represent a dressed propagator and the wavy 
lines correspond to G-matrix interactions. The first term is the Brueckner-Hartree-Fock potential while the others represent 
the 2plh/2hlp or higher contributions that are approximated through the Faddeev TDA/RPA equations. 



the above issues for 16 0. Sec. |TT] describes the formalism employed to describe the coupling of 
particle- hole (ph), particle-particle (pp) and hole-hole (hh) collective excitation to the sp propagator. 
This is done by embedding the technique of Faddeev equations into the self-consistent Green's 
function formalism 3 3|- The results of this calculation are given in Sec. IIII| where the effects of 
nuclear fragmentations and the random phase approximation (RPA) in the description of phonons are 
discussed. These calculations show that a proper description of the experimental spectral strength 
requires an improvement of the spectrum, that goes beyond the RPA approach. A first step in this 
direction is reported in Sec. IIVI where the Bethe-Salpeter equation (BSE) is employed to pursue the 
improvement of the ph-RPA approach. 



II. FADDEEV APPROACH THE CALCULATION OF THE SINGLE-PARTICLE GREEN'S FUNCTION 

We consider the calculation of the sp Green's function 



g*p{u) = — — p + ^ 



from which both the one-hole and one-particle spectral functions, for the removal and addition 
of a nucleon, can be extracted. In Eq. (JIJ), X£ = (#£ +1 |Cal*o) (3^ = (^fe _1 |c a |^)) are the 
spectroscopic amplitudes for the excited states of a system with A + l (A — 1) particles and the poles 
= E^ +l — Eq (s^ = Eq—E^~ 1 ) correspond to the excitation energies with respect to the A-body 
ground state. The one-body Green's function can be computed by solving the Dyson equation 

9 a p(u) = 9%(u) + Efl! 7 HS; s (w)^(u) , (2) 

where the irreducible self-energy T,* s (u) acts as an effective, energy- dependent, potential. The latter 
can be expanded in a Feynman-Dyson series [3, [13] in terms the exact propagator g a p(u), which 
itself is a solution of Eq. (j2J). In this expansion, Ti* s (u>) can be represented as the sum of a one- 
body Hartree-Fock potential and terms that describe the coupling between the sp motion and more 
complex excitations 0] . This separation is depicted by the diagrams of Fig. Q where the last two 
diagrams are expressed in terms of the three-line irreducible propagator R(u) which describes the 
propagation of at least 2hlp or 2plh at the same time. It is at the level of R(u) that the correlations 
involving interactions between different collective modes have to be included. 

The SCGF approach consist of first solving the self-energy and the Dyson Eq. (J2J in terms of 
a mean- field IPM propagator g^Joj). The (dressed) solution g a p{uj) is then used to evaluate the 
an improved self-energy, which then contains the effects of fragmentation. The whole procedure is 



FIG. 2: Example of diagrams that are summed to all orders by 
means of the Faddeev equations. 




iterated until self-consistency is reached. Baym and Kadanoff showed that a self-consistent solution 
of the above equation guarantees the fulfillment of the principal conservation laws 



A. Faddeev approach to the self-energy 



In the following we are interested in describing the coupling of sp motion to ph and pp(hh) 
collective excitations of the system. All the relevant information regarding the latter are included 
in the Lehmann representations of the polarization propagator 
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that describes the excited states in the system with A particles, and of the two-particle propagator 
that is relevant for the A ± 2 excitations 
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In the calculation of Sec. IIIIl these quantities have been evaluated by solving the dressed TDA/RPA 



(DTDA/DRPA) equations |20|, |2l| for the ph and pp case, respectively. The fact that a dressed 
sp propagator is used brings in to the calculation the effects of the strength distribution for particle 
and hole fragments. One must notice that the propagators Q and are in fact the solution of 
more general Bethe-Salpeter equations, that also include correlations beyond the RPA level. One of 
these is the mixing of many-phonons in the kernel of the BSE. A first calculation in this direction 
has been done for the ph propagator IT(u;) [3| and is reported in Sec. IIVI 

Once one has obtained the ph (jHJ) and pp(hh) (J3J) propagators, their coupling to the sp motion has 
to be accounted for in the nuclear self-energy. To do this the 2plh and 2hlp propagators of Fig. Q 
are computed by solving a set of Faddeev equations 0, 23[. The details of the formalism employed 
in this calculation have been discussed in Ref. UM and will not be shown here. 



For the present 

discussion it is sufficient to note that this approach computes the motion of three-quasiparticle 
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Shell TDA RPA 1st itr. 2nd itr. 3rd itr. 4th itr. 

Z Pl/2 0.775 0.745 0.775 0.777 0.774 0.776 

Z P3/2 0.766 0.725 0.725 0.727 0.722 0.724 

0.015 0.027 0.026 0.026 



Total occ. 14.56 14.56 14.56 14.57 14.58 14.63 



TABLE I: Hole spectroscopic factors (Z a ) for knockout of a £ = 1 proton from 16 0. These results refer to the initial IPM 
calculation and the first four iterations of the DRPA equations. All the values are given as a fraction of the corresponding 
IPM value. Also included is the total number of nucleons, as deduced from the complete one-hole spectral function, for each 
iteration. 



excitations in the same way it is normally done to solve the three-body problem. This allows to 
include the effects of ph and pp(hh) motion not only at the TDA level but also with more collective 
RPA phonons and beyond. Such excitations are coupled to each other by the Faddeev equations, 
generating an infinite series of diagrams like the one shown in Fig. |21 Moreover, Pauli correlations 
are properly taken into account at the 2plh/2hlp level. 



III. RESULTS FOR THE SINGLE-PARTICLE SPECTRAL FUNCTION OF le O 



In the calculations described below, the Dyson equation was solved within a model space consisting 
of harmonic oscillator sp states. An oscillator parameter b = 1.76 fm was chosen (corresponding to 
hu = 13.4 MeV) and all the first four major shells (from Is to 2plf) plus the lgg/2 where included. 
Inside this model space, the interaction used was a Brueckner G-matrix |24[ derived from the Bonn-C 
potential [Hf. The results of Refs. (HEll; suggest that this model space is large enough to properly 
account for the low-energy collective states in which we are mainly interested in, the SRC being 
accounted for through the G-matrix effective interaction. We note that energy dependence of the 
G-matrix also contributes to the normalization of the spectral strength. In Ref. it was shown 
that this accounts for the extra depletion of spectroscopic factors that is induced by SRC, at least 
for the normally occupied shells in the IPM. 

We mention that the solution of RPA equations may lead to collective instabilities whenever the 
iteration becomes too attractive. A particular sensitivity to the strength of the G-matrix interaction 
was found in the calculation of 12]| , only for the lowest isoscalar + solution. For this reason it was 
chosen to artificially shift its eigenvalue to the experimental energy of first excited state, 6.05 MeV, 
during the successive iterations. Note that the screening due the redistribution of the sp strength 
tends to reduce the instability problem in DRPA. Moreover, the origin of the instability for this 
particular solution has been identified and corrected later, in Ref. 



A. Effects of RPA correlations and fragmentation 



For an IPM ansatz the TDA calculation is equivalent to the one of Ref. M and yields the same 
result for the spectroscopic factors of the main particle and hole shells close to the Fermi energy. 
These correspond to 0.775 for pi/ 2 and 0.766 for p 3 / 2 and are reported in Table |TJ The introduction 
of RPA correlations reduces these values and brings them down to 0.745 and 0.725, respectively. 
This result reduces the discrepancy with the experiment by about 4% and shows that collectivity 
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FIG. 3: One-proton removal strength as a function of the hole sp energy ejT = Eq — -Ejf -1 for ls O for angular momentum 1 = 1 
(left) and £ — 0, 2 (right). For the positive parity states, the solid bars correspond to results for d 5 / 2 and d 3 / 2 orbitals, while 
the thick lines refer to Si/j. The top panels show the experimental values taken from 0. The mid panels give the theoretical 
results for the self-consitent spectral function. The bottom panels show the results obtained by repeating the 3 rd iteration with 
a modified ph-DRPA spectrum, in which the lowest eigenstates have been shifted to the corresponding experimental values. 



beyond the TDA level is relevant to explain the quenching of spectroscopic factors. Since the present 
formalism does not account for center-of-mass effects, the above quantities need to be increased by 
about 7% before they are compared with the experiment 26]. 

The RPA results were iterated a few times, with the aim of studying the effects of fragmentation 
on the RPA phonons and, subsequently, on the spectral strength. Since we are looking at low- 
energy excitations, it is sufficient to keep track only of the most occupied fragments that appear 
— close to the Fermi energy — in the (dressed) sp propagator, Eq. (0), while the residual strength is 
collected in a effective pole [l2l . 1271] . As shown in Table HJ only a few iterations are required to reach 
convergence. The converged distribution of one-hole strength is compared in with the experimental 
one in Fig. 01 The main difference between these results and the one obtained by using an IPM 
input is the appearance of a second smaller p 3 / 2 fragment at -26.3 MeV. This peak arises in the first 
two iterations and appears to become stable in the last one, with a spectroscopic factor of 2.6%. 
This can be interpreted as a peak that describes the fragments seen experimentally at slightly higher 
missing energy. This is the first time that such a fragment is obtained in calculations of the spectral 
strength. Further insight into the appearance of this strength is discussed in Sec. IIII Bl 

A second effect of including fragmentation in the construction of the RPA phonons is to increase 
the strength of the main hole peaks. The pi/2 strength increases from the 0.745 obtained with IPM 
input to 0.776, essentially canceling the improvement gained by the introduction of RPA correlations 
over the TDA ones. The main peak of the remains at 0.722 but the appearance of the secondary 
fragment slightly increases the overall strength at low energy as well. This behavior is due to the 
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fact that redistribution of the strength (and principally, the reduction of the main spectroscopic 
factors) tends to screen the nuclear interaction, with respect to the IPM case. Obviously, this 
makes the disagreement with experiments a little worse and additional work is needed to resolve 
the disagreement with the data. Nevertheless it is clear that fragmentation is a relevant feature of 
nuclear systems and that it has to be properly taken into account. 

Together with the main fragments, the Dyson equation produces also a large number of solutions 
with small spectroscopic factors. This strength extends down to about -130 MeV for the one-hole 
spectral function and up to about 100 MeV for the one-particle case. This background represents the 
strength that is removed from the main peaks and shifted up to medium missing energies. Another 
10% of strength is moved to very high energies due to SRC 9]. Its location cannot be explicitly 
calculated in the present approach but the effects on the reduction of spectroscopic factors at low 
energy is accounted through the energy dependence of the G-matrix. Accordingly the total number 
of nucleons deduced from the fully self-consistent spectral function is about a 10% lower than A = 16, 
as reported in Table E 

We note that the present approach gives predictions also for the strength distribution of quasipar- 
ticle excitations. Although little experimental information is available for the addition of a nucleon 
(as well as for the removal of a neutron). Recently, knockout and nucleon transfer reactions, em- 
ploying nuclear probes, have been applied to measure the hole spectroscopic factors of both stable 



and exotic nuclei |2s,|29j]. These tools can now reach better accuracies than in the past and may, in 



principle, be used to extract information on the one-particle spectral distribution. 



B. Role of and 3 excited states in O 



A deeper insight into the mechanisms that generate the fragmentation pattern can be gained by 
investigating directly the connection between the spectral function and some specific collective states. 
To clarify this point we repeated the third iteration using exactly the same input but without keeping 
the troublesome ph + state at its experimental energy, it was discarded instead. The resulting p 
hole spectral function is shown in the lower-left panel of Fig. El In this calculation no breaking of the 
p 3 /2 strength is obtained but only a single peak is found with a spectroscopic factor equal to 0.75, 
which is the sum of the two fragments that are obtained when the O2" state is taken into account. 
This result can be interpreted by considering the p%/2 fragments as a hole in the ground state and 
in the first excited + state of the 16 core, respectively. The correct fragmentation pattern can be 
reproduced only when latter two levels are close enough, in energy, to one another. 

The other two low-lying states of 16 that may be of some relevance are the isoscalar 1~ and 3~. 
These excitations are reproduced reasonably well by RPA type calculations |30| but are typically 
found at ~3 MeV above the experimental results. The lower panels of Fig. |3] show the results for the 
even parity spectral functions that are obtained when both the 3~ and the 1~ ph-DRPA solutions 
are shifted to match their experimental values. In this case, a ^5/2 hole peak is obtained at a missing 
energy -17.7 MeV, in agreement with experiment, and with a spectroscopic strength of 0.5%. 

It should be noted that many weak fragments (i.e. relative to the addition or knockout for a 
nucleon with a small spectroscopic factor) are seen throughout the nuclear table of isotopes. Most 
of these have no direct explanation in terms of mean-field sp states and contain more complicated 
components, involving 2plh/2hlp configurations and beyond. Jennings and Escher have shown that 
different definitions of the overlap wave function can be given that completely describe the sp motion 
inside the nuclear medium. Although, each of them sample the Pauli correlations in a different way, 
resulting in very different normalizations of the amplitudes of weak states This may have 

important consequences for the reliability of spectroscopic factors extracted from proton emission 
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FIG. 4: Examples of contributions involving the coupling of two 
independent ph phonons. In total, there are sixteen possible 
diagrams of this type, obtained by considering all the possible 
couplings to a ph state. The two-phonon ERPA equations sum 
all of these contributions in terms of dressed sp propagators. 
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and pick-up/stripping reactions [3 



IV. EXTENSION OF RPA AND SPECTRUM OF 16 O 

The above results suggest that the main impediment for further improvements of the description of 
the experimental spectral strength is associated with the deficiencies of the RPA (DRPA) description 
of the excited states. One important problem is the appearance of at most one collective phonon for 
a given J*", T combination while several low- lying isoscalar + and 2 + excited states are observed at 
low energy in 16 0, as well as additional 3~ and 1~ states. A possible way to proceed would be to first 
concentrate on an improved description of the collective phonons by extending the RPA to explicitly 
include the coupling to two-particle— two-hole (2p2h) states. Such an extended RPA procedure has 
been applied in heavier nuclei with considerable success Q. In order to be relevant for 16 0, this 
approach requires an extension in which the coherence of the 2p2h states is included in the form of 
the presence of two-phonon excitations (33[. This correspond to adding two-phonon configurations 
to the kernel of the BSE for the ph propagator (JHJ), as shown in Fig. |U ll3j. 

We note that there has been a tremendous progress in recent years in the microscopic de scrip tion 



scription 



of p shell nuclei using Green's Function Monte Carlo and no-core shell model methods 
A possible application of the no-core shell model to 16 O would properly include the 4p4h effects 
relevant for describing the spectrum of 16 j4£j. However the description of spectroscopic factors 
will require the construction of effective operators to include the effects of short-range correlations on 
these quantities whereas they are automatically included in the SCGF method. Whence the interest 
in studying the spectrum of 16 with the present method. 

The formalism to include two-phonon configurations has been presented in Ref. |l3j, where it is 
referred as two-phonon extended RPA (ERPA). In this work the ph-DRPA equation has been solved 
first, using the self-consistent sp propagator derived in Sec. 11111 The lowest DRPA solutions for both 
the + , 3" and 1~ channels were shifted down to their relative experimental energies and then they 
were employed to generate the two-phonon contributions for the ERPA calculation. The spectrum 
obtained for 16 is displayed in Fig. Eland Table HP together the total ph strength Z Un of each state, 

Zn, = E |<*ili<*l*tf>f • (5) 

a/3 

Table HTl also reposts the relative occupation of ph and two-phonon admixtures in the wave function 
for each solution of the ERPA equation. 

The results yield an isoscalar + state with a predominant ph character at ~17 MeV, as it was 
found in DRPA. Although this solution is now characterized by a partial mixing with two-phonon 
configurations. Table |H] shows that this is the result of mixing with the lowest solution in the 
same channel, which ends up at ~11 MeV. The latter is predominately a two-phonon state. It is 
also seen that in both cases the relevant configuration comes from the coupling of two + phonons 



8 



1- 



2~ 



\ 
1 - 



0+ 



0+ 

2- 



_L 

2- 3- iL^O* FIG. 5: Results for the DRPA and the two-phonon 



0+ — q: — :::::: "" : "2V^2+ ERPA propagator of O with the dressed input 

^ 2+ — ^" ..-•*'* 1- "" P" o+ propagator computed in Sec. Illll middle and last 

0- / —_ column respectively. In solving the ERPA equa- 

2+ / tion, the lowest 3~, 1~ and + levels of the DRPA 

— j2 — / . — 3- propagator where shifted to their experimental en- 

2: / / ergies. All other DRPA solutions were left un- 

/ / changed. The excited states indicated by dashed 

lines are those for which the (E)RPA equation 
predicts a total spectral strength Z n ^ lower than 



3- , ■' 10%. The first column reports the experimental 



0+ results |2J. 



0+ 0+ 0+ 



experiment Dressed RPA Dressed ERPA 



themselves. The calculations of Ref. [35j have shown that the bulk of the 4p4h contributions to the 
first excited state of 16 O may come from the coupling of four different phonons with negative parity 
(3~ and However, the self-consistent role of coupling positive parity states was not considered 
in that work. Both this effect, the RPA correlations and the inclusion of the nuclear fragmentation 
allow for the -at least partial- inclusion of configurations beyond the 2p2h case even when no more 
than two-phonon coupling is considered, as in this paper. A study of the importance of three- and 
four-phonon configurations will be addressed in the future. 

The low-lying negative parity isoscalar states (3~ and 1~) are only slightly affected by two-phonon 
contributions and remain substantially above the experimental energy at 9.23 and 10.90 MeV, re- 
spectively, and more correlations will be needed in order to lower their energy. We note that these 
wave functions contain a two-phonon admixture obtained by coupling the low-lying + excitation to 
either the 3~ or 1" phonons. According to the above interpretation of the first + excited state this 
corresponds to the inclusion of 3p3h and beyond. At the same time two additional negative parity 
solutions, with two-phonon character, are found at higher energy in accordance with experiment. 

The two-phonon ERPA approach also generates a triplet of states at about 12 MeV, with quantum 
numbers + , 2 + and 4 + . A similar triplet is found experimentally at 12.05, 11.52 and 11.10 MeV, 
which also corresponds to twice the experimental energy of the first 3~ phonon. The ERPA solutions 
for this triplet are almost exclusively obtained by 3~ ® 3~ states and therefore have a 2p2h character, 
in accordance with Refs. |36l |3?| • Further interaction in the pp and hh channels will be needed in 
order to reproduce the correct experimental splitting and experimental strength (35[. 



V. CONCLUSIONS 



In the present paper the Faddeev technique, combined with self-consistent Green's Function theory, 
has been applied to study 2hlp correlations at small missing energies for the nucleus of 16 0. The 
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TABLE II: Excitation energy and total spectral strengths obtained for the principal solutions of DRPA and two-phonon ERPA 
equations. The total contribution of ph and two-phonon states of the ERPA solutions are shown. For states below 15 MeV, 
the columns on the right side give the individual contributions of all the relevant two-phonon contributions. 

application of the Faddeev method allows for the first time the treatment of the coupling of ph and 
pp(hh) collective modes within a RPA framework, while the effects of fragmentation are included 
through the dressing of the sp propagator. The present approach allows to further identify the 
important role played by the low- lying excited states of 16 O. These are seen to be essential to 
generate many of the fragments with small spectrocopic factors that are seen experimentally, like 
the S1/2 and p 3/2 states of 15 N at 5.20 MeV and ~9 MeV. . 

The main impediment in obtaining a good theoretical description of the single particle spectral 
function of 16 has been identified in the poor description of the excitation spectrum, as obtained 
by solving the standard (D)RPA equations. We have improved on this by computing the effects of 
mixing of ph states with two-phonon configurations in the kernel of the Bethe-Salpeter equations. 
Moreover, configurations beyond the 2p2h level were also partially included by means of improved 
RPA and self-consistency effects. The results show that these contributions explain the formation 
of several excited states observed at low energy which are not obtained by RPA calculations. Still it 
appears that a full solution of the spectrum of 16 with this method requires to consider the mixing 
of up to four phonons and the interaction in the pp and hh channels |35|. The inclusion of such 
effects will be the topic of future work. 
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